Theory of Magnetic Anisotropy in IIIi.^Mn^V Ferromagnets 
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We present a theory of magnetic anisotropy in IIIi-ajMn^V diluted magnetic semiconductors with 
carrier-induced ferromagnetism. The theory is based on four and six band envelope functions models 
for the valence band holes and a mean-field treatment of their exchange interactions with Mn ++ 
ions. We find that easy-axis reorientations can occur as a function of temperature, carrier density p, 
and strain. The magnetic anisotropy in strain-free samples is predicted to have a p 5 ^ 3 hole-density 
dependence at small p, a dependence at large p, and remarkably large values at intermediate 
densities. An explicit expression, valid at small p, is given for the uniaxial contribution to the 
magnetic anisotropy due to unrelaxed epitaxial growth lattice-matching strains. Results of our 
numerical simulations are in agreement with magnetic anisotropy measurements on samples with 
both compressive and tensile strains. We predict that decreasing the hole density in current samples 
will lower the ferromagnetic transition temperature, but will increase the magnetic anisotropy energy 
and the coercivity. 
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I. INTRODUCTION 

The discovery of carrier-mediated ferromagnetismErll 
in IIIi-^Mn^V and dopedcl IIi-^Mn^VI diluted mag- 
netic semiconductors (DMS's) has opened up a broad 
and relatively unexplored frontier for both basic and 
applied research. ExperimentsBB in Gai_ x Mn x As and 
Ini-xMn^As have demonstrated that these ferromagnets 
have remarkably square hysterisis loops with coercivities 
typically ~ 40Oe, and that the magnetic easy axis is de- 
pendent on epitaxial growth lattice-matching strains. In 
this paper we discuss the magnetic anisotropy proper- 
ties of IIIi-^Mn^ DMS ferromagnets, predicted by a 
mean-field-theoryB of the exchange interaction coupling 
between localized magnetic ions and valence band free 
carriers. We use phenomenological four or six band en- 
velope function models, depending on the carrier density 
p, in which the valence band holes are characterized by 
Luttinger , spin-orbit splitting, and strain-energy param- 
eters. 

The physical origin of the anisotropy energy in our 
model is spin-orbit coupling in the valence band. Our 
work is based in part on theoretical descriptions devel- 
oped by GajEI and BastardB to explain the optical prop- 
erties of undoped, paramagnetic DMS's. As the criti- 
cal temperature is approached, the mean-field-theory we 
employO reduces to an earlier theoryu that invokes gener- 
alized RKKY carrier-mediated interactions between lo- 
calized spins. The two approaches differ, however, in 
their description of the magnetically ordered state. As 
this work wasj-nearing completion we learned of a closely 
related studyH which uses the same mean-field theory 
to address critical temperature trends in this material 



class and which also comments on the mean-field the- 
ory's ability to address magnetic anisotropy physics. We 
are aware of three elements of the physics of these mate- 
rials which make the predictions of our mean-field theory 
uncertain: i) we do not account for the substantial dis- 
order which is usually present in these ferromagnets; ii) 
we do not include effects due to interactions among the 
itinerant holes and iii) we do not account for correla- 
tions between localized spin-configurations and itinerant 
hole states. The importance of each of these deficien- 
cies is difficult to judge in general, and probably depends 
on adjustable material parameters. In our view, it is 
likely that there is a substantial range in the parameter 
space of interest where the predictions of the present the- 
ory are useful. We expect that important progress can 
be made by comparing this simplest possible theory of 
carrier-induced DMS ferromagnetism with experiment. 

This work has two objectives. Most importantly, we 
have attempted to shed light on how various adjustable 
material parameters can influence magnetic anisotropy. 
Secondarily, we have made an effort to estimate the mag- 
netic anisotropy energy in those cases where experimen- 
tal information is presently available. Our hope here is 
to initiate a process of careful and quantitative compari- 
son between mean- field theory and experiment, partially 
to help judge the efficiency of this approximation in pre- 
dicting other physical properties. Even in the mean-field 
theory, we find that the magnetic anisotropy physics of 
these materials is rich. We predict easy axis reorienta- 
tions as a function of hole density, exchange interaction 
strength, temperature, and strain and identify situations 
under which IIIi-^Mn^V ferromagnets are remarkably 
hard. 
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In Section II we detail our mean-field-theory of the 
ordered state. The theory simplifies in the limit of low- 
temperature and low hole densities. Our results for this 
limit, presented in Section III, predict a (111) easy axis 
in the absence of strain, and a magnetic anisotropy en- 
ergy which is approximately 10% of the free carrier band- 
energy density. This value is extremely large for a cubic 
metallic ferromagnet; typical ratios in transition metal 
ferromagnets are smaller than 1CP 6 for example. The 
anisotropy energy in this limit varies as the free-carrier 
density to the 5/3 power and is independent of the ex- 
change coupling strength. Explicit results for the strain- 
dependence of the magnetic anisotropy in the same limit 
are presented in Sec. IV. We find that unrelaxed lattice- 
matching strains due to epitaxial growth contribute a 
uniaxial anisotropy which favors magnetization orienta- 
tion along the growth direction when the substrate lattice 
constant is larger than the ferromagnetic semiconductor 
lattice constant and an in-plane orientation in the oppo- 
site case. Unfortunately, perhaps, the simple low density 
limit does not normally apply in situations where high 
critical temperatures are expected. The more compli- 
cated, and more widely relevant, general case is discussed 
in Section V. We find that magnetic anisotropy has a non- 
trivial dependence on both temperature and exchange 
coupling strength and that easy axis reversals occur, in 
general, as a function of either parameter. According 
to our theory, anisotropy energy densities comparable to 
those in typical metallic ferromagnets are possible when 
the exchange coupling is strong enough to depopulate all 
but one of the spin-split valence bands, even with satura- 
tion magnetization values smaller by more than an order 
of magnitude. In the limit of large hole densities, we 
find that the anisotropy energy of strain-free samples is 
proportional to hole density p^ 1 and exchange coupling 
Jp d . We find that in typical situations a strain eo of only 
~ 1% is sufficient to overwhelm the cubic anisotropy of 
strain-free samples. We conclude in Section VI with a 
discussion of the implications of these calculations for 
the interpretation of present experiments, and with some 
suggestions for future experiments which could further 
test the appropriateness of the model used here. 

II. FORMAL THEORY 

Our theory is based on an envelope function descrip- 
tion of the valence band electrons, and a-scin representa- 
tion for their, kinetic- exchange inteactionEj with localized 
d electronsEil on the Mn ++ ions: 

n = n m + n b + Jpd s(n - Ri), (i) 

i,I 

where i labels a valence band hole and / labels a mag- 
netic ion. In Eq. ([!]) , TL m describes the coupling of mag- 
netic ions with total spin quantum number J = 5/2 



to an external field (if one is present), Si is a local- 
ized spin, Si is a hole spin, and Jib is either a four or 
a six-band envelope-function HamiltonianEa for the va- 
lence bands. In this paper we do not consider exter- 
nal magnetic fields so H m — > 0. The four-band Kohn- 
Luttinger model describes only the total angular momen- 
tum j = 3/2 bands, and is adequate when spin-orbit cou- 
pling is large and the hole density p is not too large. As 
discussed later, in the case of GaAs, a four-band model 
suffices for p < 10 18 cm" 3 . In IIIi-^Mn^V semiconduc- 
tors, the four j = 3/2 bands are separated by a spin-orbit 
splitting A so from the two j — 1/2 bands. In the rele- 
vant range of hole and Mn ++ densities, no more than 
four bands are ever occupied. Nevertheless, mixing be- 
tween j = 3/2 and j — 1/2 bands does occur, and it 
can alter the balance of delicate cancellations which often 
controls the net anisotropy energy. The exchange inter- 
action between valence band holes aiid localized moments 
is believed to be antiferromagBstwEJ, i.e. J p d > 0. For 
GaAs, experimental estimatestLiatj of J p d fall between 
0.04eVnm 3 and 0.15eVnm 3 , with more recent work sug- 
gesting a value toward the lower end of this range. 

The form of the valence band for Bloch wavevec- 
tors near the zone center in a cubic semiconductor 
follows from k_i. p perturbation theory and symmetry 
considerations. til The four band (j — 3/2) and six band 
(j = 3/2 and 1/2) models are known as Kohn-Luttinger 
Hamiltonians and their explicit form is given in the Ap- 
pendix. The eigenenergies are measured down from the 
top of the valence band, i.e., they are hole energies. 
The Kohn-Luttinger Hamiltonian contains the spin-orbit 
splitting parameter A so and three other phenomenolog- 
ical parameters, 71, 72 and 73. These are accurately 
known for common semiconductors. For GaAs and In As, 
the two materials in which IIIi_ x Mn :E V ferromagnetism 
has been observed, A so = 0.34 eV and 0.43 eV, and 
(71,72,73) = (6.85,2.1,2.9) and (19.67,8.37,9.29) re- 
spectively. Most of the specific illustrative calculations 
discussed below are performed with GaAs parameters. 

Our calculations are based on the Kohn-Luttinger 
Hamiltonian and on a mean-field theory in which corre- 
lations between the local-moment configuration and the 
itinerant carrier system are neglected. We comment later 
on limits of validity of this approximation. There are 
a number of equivalent ways of developing this mean- 
field theory formally. In the following paragraphs we 
present a view which is convenient for discussing mag- 
netic anisotropy. 

In the absence of an external magnetic field, the par- 
tition function of our model may be expressed exactly as 
a weighted sum over magnetic impurity configurations 
specified by a localized spin quantization axis, M, and 
azimuthal spin quantum numbers to,/: 

Z = ]Texp(-F b [TO/]/fc/jT), (2) 

JTi I 

where Fb[m{\ is the valence band free energy for holes 
which experience an effective Zeeman magnetic field 
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h(r)[mi] = -JpdMy^mi 6(r- Ri 



(3) 



The mean- field theory consists of replacing h(r)[mj] by 
its spatial average for each magnetic impurity config- 
uration, thereby neglecting correlations between spin- 
distributions in local-moment and hole subsystems. The 
effective Zeeman magnetic field experienced by the holes 
then depends only on M, the direction of the local- 
moment orientation, and the mean azimuthal quantum 
number averaged over all local moments, M: 



h MF (M) = J pd N Mn MM = hM, 



(4) 



where Nmu = Ni/V is the number of magnetic impuri- 
ties per unit volume. The mean-field partition function 
is 



Z MF (M) = exp((7V / Ts(M) - F b {h))/k B T), 
where the entropy per impurity is defined by 



s(M) 



k B lim 



(5) 



(6) 



and F b (h) is the free-energy of a system of non-interacting 
fermions with single-particle Hamiltonian H b — hM ■ s. 

Following standard 'large number' arguments s(M) is 
readily evaluated by considering an auxiliary system con- 
sisting of magnetic impurities coupled only to an exter- 
nal magnetic field H. For this model problem, a familiar 
exerciseEHl gives the result 



M{H) = JBj(x) 



(7) 



where x = gL^ B HJ/k B T, g^ is the Lande g-factor of the 
ion, fis is the electron Bohr magneton, and 

Bj(x) = 2J ^ 1 coth [(2J + l)x/2J] - -^j coth(x/2J), 

(8) 



is the Brillouin function. The Brillouin function is a one- 
to-one mapping between reduced fields x in the interval 
[0, oo] and reduced magnetizations M / J in the interval 
[0,1]; the inverse function BJ 1 maps M/J to x. Since 
the magnetization maximizes s{M) + gL^BHM/ksT, 



ds(M) 
dM 



-g L li B H/k B T. 



(9) 



Eq. (|^) can be used to eliminate H and arrive at the 
following explicit expression for s(M): 



s(M) = ki 



dx x 



Bj L (M/J) 



dBjjx) 
dx 



(10) 



The J = 5/2 result for s(M) is illustrated in Fig. |l[ The 
entropy per impurity vanishes for M = J = 5/2 because 
there is a single configuration with X)/ TO ^ = NjJ, and 
approaches In (2 J + 1) w 1.79 for M -> 0. 
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FIG. 1. Entropy per localized spin as a function of average 
polarization. 

The mean-polarization of the localized spins at a given 
temperature and for a given orientation of the local mo- 
ments is determined by minimizing the mean-field free 
energy 

F MF (M) = -k B T\nZ MF {M) 

= F b (h = NunJpdMM) - k B TNis(M), (11) 

with respect to M. Setting the derivative to zero gives 



ds(M) _ J pd dF b (hM) 
dM ~ k B TV dh 



(12) 



Comparing with Eq. (g), it follows that Fmf(M) is min- 
imized by M = JBj(g L fi B H e ffJ/k B T) = JBj(x eff ) 
where 



X eff 



gL^sHeff _ J pd dF b {hM) 
k B T ~ Vk B T dh ' 



(13) 



It follows that h = NMnJpdM is determined by solving 
the self-consistent equation 



h = N Mn J P dJBj [x e ff{h)\ . 



Note that 



dF b (hM) 
dh 



= (S tot -M), 



(14) 



(15) 



where Stot is the total hole spin and the angle brackets 
indicate a thermal average for the non-interacting valence 
band system. Since the valence band system experiences 
an effective Zeeman coupling with strength proportional 
to h and direction —M, it is clear that the right hand 
side of Eq. ( |l5| ) is negative in sign and that its magnitude 
increases monotonically with h, making it easy to solve 
Eq. (114) numerically. 
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To simplify the calculations presented in subsequent 
sections, we take advantage of the fact that temperatures 
of interest are almost always considerably smaller than 
the itinerant carrier Fermi energy. This allows us to re- 
place Fb(h) by the ground state energy Eb(h). Then, us- 



ing Eq. (y_3J) and Eq. (Jl l[) , a single calculation of Eb(hM) 
over the range from h = to h = NMnJpdJ may be used 
to determine the local- moment magnetization M(T) and 
the free energy F(T) = Fmf{M{T)) at all temperatures. 

The mean-field theory critical temperature can be 
identified by linearizing the self-consistent equation at 
small h. We find that 



k B T c (M) = 



J {J + 1) N Mn J 2 d d 2 F b (hM) , 



V 



dh 2 



h=0- 



(16) 



The second derivative of the valence band free-energy 
with respect to field is proportional to its Pauli spin- 
susceptibility, which is in turn proportional to the va- 
lence band density of states at the Fermi energy, and to 
p 1 / 3 at small p. In the absence of strain, it follows from 
cubic symmetry that the right-hand-side of Eq. (^) is 
independent of M. Below the critical temperature how- 
ever, the mean-field free energy does depend on M; this 
dependence is the magnetic anisotropy energy we wish 
to calculate. We will see that the dependence of the 
anisotropy energy on hole density is very different from 
that of the critical temperature. 



III. MAGNETIC ANISOTROPY IN THE 
STRONG EXCHANGE COUPLING LIMIT 

Our mean-field theory simplifies at low temperatures 
and, for the four-band model, simplifies further when h is 
much larger than the characteristic energy scale of occu- 
pied Kohn-Luttinger states. A convenient typical energy 
scale is the h = hole Fermi energy e fo ■ For a given value 
of NMnJpd the largest value of h is reached at T = 0. 
Then, since H e ff is always non-zero, x e ff — » oo and the 
solution to the mean-field equations is M — J, implying 
that h = NMnJpd J for every orientation M, At T = 
there is no entropic contribution to the free-energy and 



F MF (T = 0) = E b (h = NMnJpd J M). 



(17) 



We note in passing that the magnetization density at 
T = has contributions from the localized spins and the 
itinerant spins: 



M S (T = 0) = 2fi B JN Mn + -^{Stot ■ M), 



(18) 



where k is an additional parameter of the Luttinger 
model. Because of the antiferromagnetic exchange in- 
teraction, the two terms here will tend to have opposite 
signs with the first term, which is independent of the 
hole density, typically very dominant. For h 3> e^o two 
further simplifications occur. When the splitting of the 



hole bands by the effective Zeeman coupling is sufficiently 
large, or the hole density p is sufficiently small, only the 
lowest energy hole band will be occupied. Furthermore, 
the effective Zeeman-term will dominate the mean-field 
single-particle Hamiltonian and, as we detail below, the 
envelope function spinor for this occupied hole state has 
a simple analytic expression. In this section we assume 
that the spin-orbit splitting energy A so is much larger 
than all other energies so that we can work with a four 
band model. More generally, the anisotropy will depend 
on h/ A so , even in the limit of small hole densities. 

To judge whether or not this limit can be achieved 
in practice, we estimate the .Fermi energy of holes using 
the spherical approximatiorJ13 in which the doubly de- 
generate h = bulk heavy hole and light hole bands are 
parabolic with masses rtih = m/(7i — 27) ~ 0.498m and 
mi = 771/(71 + 27) ~ 0.086m respectively. An elementary 
calculation then gives 



h 2 ,37T 2 /K2/3 



7o, 



2m y 2 

where n = N^/V is the free-carrier density and 



7o 



.( 7 i-27)- 3 / 2 + (7i+27)- 3/2 - 



-2/3 



(19) 



(20) 



For GaAs 70 ~ 3.05: the Fermi energy is the same as that 
of a system with four identical effective mass m/70 bands. 
Typical high T c Gai-zMn^As ferromagnetic semiconduc- 
tor samples have p ~ O.lnm -3 and Nmh ~ l.Onm -3 
(a; ~ 0.05), although these parameters can presumably 
be varied widely. Choosing a J p d value in the mid-range 
of estimated values (~ 0.006Ry ran 3 ) these parameters 
imply that h ~ 0.015Ry and e^o ~ O.OlRy. h is nei- 
ther large compared to e^o, nor small compared to A so . 
Thus, the simple expressions discussed in this section are 
not accurate for current high T c systems. As our abil- 
ity to engineer materials improves it should, however, be 
possible to grow samples which are in the limit discussed 
here. Since h is comparable to e^o, we know, even be- 
fore performing detailed calculations, that valence band 
quasiparticle spectra in paramagnetic and ferromagnetic 
states will differ qualitatively. 

In the large h limit the lone occupied spinor at each 
wavevector k will be the member of the j = 3/2 mani- 
fold for which the total angular momentum is aligned in 
the direction — M; the origin of the minus sign here is 
the antiferromagnetic nature of the interaction between 
localized spins and hole spins. Explicit expressions for 
the expansion of such a spin coherent states in terms of 
the eigenstates of j z are knownEiJ: 



M) = u 3 |3/2) + V3[u 2 w|l/2) + uv 2 \ 
+ v 3 \ - 3/2 > . 



1/2)] 



(21) 



In Eq. (Elj) u = -iexp(-«0/2)sin(#/2) and r 
2)co 



and 

i exp(i<j)/2~)cos(8/2) where 9 and <f> are the spherical co- 
ordinates which specify the unit vector —M. In the large 
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h/eFo limit, the band term in the single-particle Hamil- 
tonian may be treated using first order perturbation the- 
ory. Taking the expectation value of the Kohn-Luttinger 
Hamiltonian in the spin coherent state we find that 



e(k) 



h 
2 
h 

2 + ~2m 



(-M\H L (k)\- 
■7(M,A). 



M) 



(22) 



The first term on the right handjiide of Eq. ( p^ ) reflects 
the spin coherent state propertyE3, M ■ s\ — M) = — \ — 
M)/2. In Eq. (^2|) we have noted that for any M and any 
fc, the dependence of hole energy on k = \k\ is quadratic. 
Using this property, it follows that the Fermi energy 



^(M) = — (6tt 2 p) 2 / 3 7(M), 



(23) 



and that the ferromagnetic ground state energy density 
is 



E b (M) 
V 



hp 



tP e F (M), 



where 



<lk 



7(M)= [ / — ( 7 (M,fc)) 



-3/2] -2/3 



(24) 



(25) 



In analogy with the h = quantity 70 defined in 
Eq. (|2(]), 7(M) is an average of the band energy curva- 
ture over reciprocal space directions fc, with the smaller 
values of j(M, k) weighted more heavily. Note that the 
factor 37r 2 /2 in Eq. ( fl9|) is replaced by 67r 2 in Eq. ( |23| ) 
because only one band is occupied in this limit, instead 
of the four which are occupied at h = 0. m/7(M) may 
be thought of as a spin-orientation dependent effective 
mass, which is readily evaluated as a function of M, given 
the Luttingcr parameters of any material. Although the 
magnetic condensation energy has a term proportional 
to J p d, only the band-energy contributes to the M de- 
pendence of the ferromagnetic ground state energy. The 
magnetic anisotropy energy in this limit is independent 
of Jpd and proportional to the hole density p to the 5/3 
power. 

We have evaluated 'y(M) as a function of angle for 
the Luttinger parameters of GaAs and InAs. As dis- 
cussed in more detail later, we always find that mag- 
netic anisotropy in the absence of strain is well described 
by a cubic harmonic expansion truncated at sixth opder, 
an approximation commonly used in the literatureEa on 
magnetic materials. The corresponding cubic harmonic 
expansion for 7(M) is 

7(M) = 7«100» + ll a (M 2 x M 2 v + M 2 M 2 + M 2 X M 2 Z ) 



+ 72 ca M 2 M 2 M 2 . 



(26) 



In our calculations we find that ~f{M) < 71 for all direc- 
tions M. This property reflects the fact that small cur- 
vature (large Fermi wavevector) directions are weighted 
more highly in calculating the total hole energy. For both 
InAs and GaAs we find that the dominant fourth order 
cubic anisotropy coefficient, 7™ < 0, indicating nickel- 
type anisotropy with easy axes along the (111) cube di- 
agonal directions. At a qualitative level, the source of 
the higher total hole energy when the moment orienta- 
tion is along a (001) (cube edge) direction is easy to un- 
derstand. With such a moment orientation, the occupied 
hole orbital has j z — —3/2 and energy dispersion given 
by H hh (k) in Eq. ©. It follows that -y(M = (100), k) 
has the relatively large value 71 + 72 for all orientations 
of k in the x — y plane. These large values of "f(M, k) are 
important in the average and cause the average over k to 
reach its maximum for this orientation of M . 

The cubic magnetic anisotropy energy coefficients in 
this limit are given by 



Kf a = 



3 hf_ 



(6tt 2 p) 2 / 3 7™. 



(27) 



Values of 7™ for GaAs and InAs are listed in Table 
I. We will see later that these expressions apply up to 
p ~ 10 18 cm~ 3 . Inserting this value for the hole density 
gives coefficients ~ 2kJm~ 3 for GaAs host material and 
~ 4kJm~ 3 for InAs host materials; magnetic anisotropy 
is twice as strong in InAs in the strain free case. These 
anisotropy energy coefficients are not so much smaller 
than those of the cubic metallic transition metal ferro- 
magnets, despite the much higher carrier densities in 
the metallic case. We will see below that the scale of 
the semiconductor magnetic anisotropy energy does not 
change as the carrier density increases from 10 18 cm -3 
to ~ 10 21 cm~ 3 . The relatively large anisotropy energies 
occur despite the fact that the saturation moments M s 
of IIIi-^Mn^V ferromagnets are more than an order of 
magnitude smaller than their cubic metal counterparts. 
It follows from these values that the magnetic hardness 
parameters of the IIIi-^Mn^V ferromagnets, 



AT 



1 1/2 



(28) 



will typically be larger than one. This is unusual in cubic 
materials and occurs because spin-orbit coupling has a 
much stronger influence on semiconductor valence bands 
than on transition metal d-bands. 

As we discuss at length in Section V, magnetic 
anisotropy does not continue to increase rapidly with 
hole density once two or more bands are occupied in 
the metallic state. In current high T c samples, we will 
find that several bands are always occupied, even at zero 
temperature. The simple limit discussed in this section 
demonstrates that anisotropy energies, T = saturation 
moments, and critical temperatures will have radically 
different dependencies on engineerablc parameters. 
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IV. STRAIN DEPENDENCE OF MAGNETIC 
ANISOTROPY: LOW HOLE DENSITY LIMIT 

Because of the low soluability of Mn in III-V semi- 
conductors, IIIi-^Mn^V materials with x large enough 
to produce cooperative magnetic effects cannot be ob- 
tained by equilibrium growth. The MBE growth tech- 
niques which have been successfully developed^ produce 
IIIi_ x Mn. E V films whose lattices are locked to those of 
their substrates. X-ray diffraction studied 3 ! have estab- 
lished that the resulting strains are not relaxed by dislo- 
cations or other defects, even for thick films. Strains in 
the IIIi-^Mn^V film break the cubic symmetry assumed 
in the previous section. Fortunately, the influence of 
MBE growth lattice-matching strains on the hole bands 
of cubic semiconductors is wellt3E-a understood. For the 
(001) growth direction used to create IIIi-^Mn^V films, 
strain generates a purely diagonal contribution to the 
four band single-particle envelope function Hamiltonian 
in the representation we use in this paper, adding contri- 
butions 5th and 5ei respectively to j z — ±3/2 heavy hole 
and j z = ±1/2 light hole entries. The energy shifts are 
related to the lattice strains byO 



= [ - 2ai(C u - C 12 ) - ^-(C n + 2Ci 2 )] , (29) 
On ^ 



a 2 , 



Sei = [ - 2 0l (Cu - Ci 2 ) + ^(Cn + 2C 12 )] , (30) 

where eo is the in-plane strain produced by the substrate- 
film lattice mismatch: 



a 2 , 



eo 



as ~ a-F 



(31) 



In Eqs. (]29|-|31|), the are the elastic constants of the 
unstrained IIIi-^Mn^V film, which we will assume to be 
identical to those of the host III-V material, as is the lat- 
tice constant of the substrate on which the IIIi-^Mn^V 
film is grown, clf is the unstrained lattice constant of 
bulk IIIi-^Mn^V, and a\ and a 2 are phenomenological 
deformation potentials whose values for common III-V 
semiconductors are known. For the six band model the 
strain Hamiltonian includes off-diagonal elements and is 
given up to a constant times the unit matrix by 



H, 



strain 



= Ten 



( 




















1 
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v° 


1 

V2 











1 

2 



(32) 



where V = e; — eh — a 2 eo(Cn+2Ci 2 )/Cn. For GaAs>and 
InAs, T = -0.2382 Ry and -0.2762 Ry respectively.^ 



As in the previous section, we can derive an explicit 
expression for the strain contribution to the magnetic 
anisotropy energy when h ^> cfo an d h ^> T, allowing 
band and strain terms to be treated as a perturbative 
correction to the effective Zeeman coupling, and h <C 
A so , allowing a four-band model to be used. In this way 
we obtain 

Estrain(M) P 3p COS 2 6 

= ~[35e h + 8ei] + [Se h -Sei\. (33) 

Strain produces a uniaxial contribution to the mag- 
netic anisotropy of IIIi-^Mn^V films, that favors ori- 
entations along the growth direction when strain shifts 
heavy holes down relative to the light holes and orienta- 
tions in the plane in the opposite circumstance. Using 
Eq. (|2^) and Eq. ^0|) and deformation potential elas- 
tic constant valuestJ^we find that a contribution to the 
energy density given by K str ain sin 2 (#) where the uniax- 
ial anisotropy constant is K stra in = — 0.36Ry eo p for 
GaAs and K stra in = — 0.41Ry eo p for InAs. Since 
a 2 < 0, compressive strain (eo < 0) lowers the heavy 
hole energy relative to the light hole energy and favors 
moment orientations in the growth direction while tensile 
strain (eo > 0) favors moment orientations perpendicu- 
lar to the growth direction. Since the lattice constant of 
Ga,i^ x Mn x As is largerB than that of GaAs, while that of 
Ini^Mn-r As is smaller than that of InAs, Gai-^Mn^As 
on GaAs is under compressive strain and Inj-jMn^As 
on InAs is under tensile strain. AssummingB Vergard's 
law, e = .0004 and e = -.0028 for InAs and GaAs 
respectively at x = 0.05. For p — 10~ 3 nm~ 3 , a density 
for which this low-density result is still reasonably reli- 
able, we find that K strain = — 0.36kJm -3 for InAs and 
K s train = 2.2kJm -3 for GaAs. At these densities, strain 
and cubic band contributions are comparable in GaAs, 
but the latter contribution is dominant in the InAs case. 

Strain-induced anisotropy energies can be modified by 
growing the magnetic film on relaxed buffer layers. This 
effect has been demonstratedEj by Ohno et al. who 
showed that the easy axis for Gai-^Mn^s films changes 
from in-plane to growth-direction when the substrate is 
changed from GaAs to relaxed (In,Ga)As buffer layers. 
For 15% In, the magnetic film strain changes from com- 
pressive eo = —0.0028 to tensile eo = 0.0077 when this 
change is made. We note that the sense of this change 
is opposite to that predicted by our large h, small p ana- 
lytic result which predicts that compressive strains favor 
growth direction orientations. Similarly, Ini-^Mn^As on 
InAs is under a small tensile strain but is observed to 
have a growth direction easy axis. (Recall that in the 
large h limit the band anisotropy makes the growth di- 
rection the hard axis.) This distinction may be taken 
as an experimental proof that several hole bands are oc- 
cupied in the magnetic ground state of these materials. 
The strain anisotropy energy does change sign at smaller 
values of h partly because the first band to be depop- 
ulated has primarily heavy-hole character. We will see 
in the next section that the mean-field does predict the 



G 



correct sign for the strain anisotropy energy at exper- 
imental hole densities. We conclude from the present 
considerations that strain can play a strong role in band- 
structurc-enginccring of IIIi-^Mn^V ferromagnet mag- 
netic properties. 



V. PARTIALLY POLARIZED HOLE BAND 
STATES: MAGNETIC ANISOTROPY IN THE 
GENERAL CASE 

The results for Fb(hM) discussed in the previous sec- 
tion become accurate when the effective Zeeman cou- 
pling h is large enough to reduce the number of occu- 
pied hole bands to one, and the Fermi energy remains 
safely smaller than the spin-orbit splitting. The situa- 
tion is much more complicated at smaller h and larger 
hole densities. Then several bands are occupied, even at 
T = 0, and these usually give competing contributions 
to the magnetic anisotropy. It is not true in general that 
band and strain contributions to the magnetic anisotropy 
are simply additive. We expect that it will eventually be 
possible to realize a broad range of material parameters, 
and hence a broad range of h values, in llli^ x Mn x V fer- 
romagnets. The range of possibilities is immense and 
accurate modeling of a particular sample will require ac- 
curate values for p, J p d, and Nmu for that material. 

In this section we discuss a series of illustrative cal- 
culations, starting with ones performed using a four- 
band model of strain-free Gai-^Mn^As at hole density 
p = 0.1nm~ 3 . The valence band energy density is 



Eb 
V 



e F - ej(k, h))ej(k, h), 



(34) 



k j=l 



where €f is the Fermi energy, k is the Bloch wavevector, 
and the mean- field-theory quasiparticle energies 6j (k, h) 
are eigenvalues of the Nf, x Nb single-particle Hamilto- 
nian (Nb is the number of bands included in the envelope 
function Hamiltonian) 



H h =H L + h-s + H s 



strain ■ 



(35) 



For a strain free model we set -Hstrain = 0. In Fig. g we 
plot the calculated spin-polarization per hole as a func- 
tion of h for a growth direction field orientation; recall 
that this quantity can be obtained by differentiating the 
energy per hole with respect to h and that the effective 
field seen by the localized spins is obtained by multiply- 
ing this quantity by J p d p. The hole spin polarization 
increases linearly at small h with a slope proportional to 
the valence band Pauli susceptibility. We see that for the 
hole density of Fig. [| complete hole spin-polarization is 
approached only at values of h comparable to or larger 
than the spin-orbit splitting of GaAs, so that the four 
band model is not physical in this regime. For any given 
model, and a given moment orientation, a single calcula- 
tion of this type provides all the microscopic information 



required to solve the mean-field equations at all temper- 
atures. 





-0.1 



■£ -0.2 

co -0.3 

v 

-0.4 
-0.5 



0.2 0.4 

h(Ry) 



0.6 



FIG. 2. Valence band spin per hole as a function of effective 
Zeeman field strength h. 

For a fixed values of J p d, Nmu, and p, h must be eval- 
uated as a function of temperature by solving the self- 
consistent mean- field equation, Eq. (|14|), and using nu- 
merical results like those plotted in Fig. g. Results for 
h(T) calculated for J p d — 0.15eVnm 3 , Nmu = limi -3 , 
and p = O.lnm -3 are illustrated in Fig. [| for high sym- 
metry moment orientation directions. For these param- 
eters the critical temperature T c (h(T) = for T > T c ) 
is ~ 100K, in rough agreement with experiment. There 
are, however, other choices of parameters which are also 
consistent with the measured critical temperature. In 
addition, as we discuss further in the next section, it is 
not clear that this level of theory should always yield 
accurate results for T c . 
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FIG. 3. The dependence of the effective Zeeman coupling 
strength seen by the magnetic impurities h on temperature 
for the high symmetry directions < 100 >, < 110 >, and 
< 111 >. These results were calculated using a strain free 
four-band model for J p d = 0.15eVnm 3 and Nm„ = l.Onm -3 . 
h(T) is determined by solving the self-consistent mean-field 
equation. 

At T = 0, h reaches its maximum value, JpdNninJ ~ 
0.0275Ry. The dependence of energy on magnetization 
orientation at this value of h is illustrated in Fig. || and 
compared with the cubic harmonic expansion truncated 
at 6th order. The coefficients of this expansion are fixed 
by energy per volume calculations in (100), (110), and 
(111) directions. In Fig. EL and in all other cases we have 
checked, the truncated cubic harmonic expansion is very 
accurate. It is therefore sufficient to evaluate the energy 
per volume in the high-symmetry directions and to use 



4(^(110) -£,,(100)) 
V 

27£ 6 (111) - 36£ b (110) + 9£ b (100) 
V ' 



(36) 



Results for Kf a (h) obtained for the four band model in 
this way are summarized in Fig. 0. These results can be 
combined with those in Fig. || to obtain the model's cu- 
bic anisotropy coefficients as a function of temperature, 
hole density and h. Here we see explicitly the anisotropy 
reversals which commonly accompany hole band depop- 
ulations. We note in Fig. || that the analytic results of 
Section III are recovered only for very large values of h 
at this density. 
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FIG. 4. Hole energy per particle in units of h as a function 
of magnetization orientation for h — 0.0275Ry). The solid line 
is a cubic harmonic expansion fit to these results truncated 
at 6th order. At this value of h, three hole bands are partly 
occupied and the easy axes are the cube edge directions. 
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FIG. 5. The dependence of the crystalline anisotropy co- 
efficients, Ki a (h) and (h) on h for a four-band model with 
p = O.lnm -3 . 

The valence band Fermi surfaces of IIIi-^Mn^V fer- 
romagnets will be strongly dependent on both tempera- 
ture and moment direction orientation. Four-band model 
Fermi surface intersections with the k z = plane are 
illustrated in Figs. ||-|| for p = O.lnm -3 at h = 0, 
h = O.OlRy, and h = 0.0275Ry respectively. All figures 
are for moments oriented in the (100) direction. The 
largest of these three values of h is the T = effective 
field (JpdNMnJ) for x = 0.05 and J p d at the high end 
of literature estimates. (J p d = 0.15eVnm 3 .) These three 
figures represent the mean-field-theory Fermi surfaces at 
three different temperatures. In the spherical approxi- 
mation, the h = Fermi energy at this hole density is 
€fo ~ O.OlRy, so a strong distortion of the bands is ex- 
pected in the magnetic state. At h = 0, the hole bands 
occur in degenerate pairs; we refer to the two less dis- 
persive bands which occupy the larger area in k space as 
heavy-hole bands, although this terminology lacks precise 
meaning in the general case. As h increases, both heavy 
and light hole bands split. For small h, the minority-spin 
heavy-hole band occupation decreases rapidly and all 
other band occupations increase. The heavy-hole minor- 
ity spin band is completely depopulated for h ~ 0.04Ry. 
Once this band is empty, the light-hole minority-spin 
Fermi radii begin to shrink rapidly. As seen in Fig. ||, 
this second band is nearly depopulated at h = 0.0275Ry. 
At still stronger fields, the majority-spin light-hole band 
is depopulated and the single-band limit addressed in 
preceding sections is achieved. For p = O.lnm" 3 , the 
single-band limit is achieved only at h values for which 
the four-band model breaks down. 
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FIG. 6. Fermi surface intersection with the k z = plane 
for a four-band model with p — O.lnm" 3 and h — 0. Doubly 
degenerate heavy-hole (larger contour) and light-hole (smaller 
contour) bands are occupied. 
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FIG. 7. Intersection of the Fermi surface and the k z = 
plane for a four-band model with p = O.lnm -3 and 
h — O.OlRy. This value of h solves the mean-field equations 
equation at T = OK, Nitn = lnm~ 3 , and J p d f=H0.05eVnm 3 , 
which is near the lower experimental estimates for the ex- 
change coupling constant. (For J p d = 0.15eVnm 3 this value 
of h solves the mean- field equations at T = 85K). The mag- 
netization orientation is along the (100) direction. Heavy-hole 
and light-hole bands are split at non-zero h. 
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FIG. 8. Intersection of the Fermi surface and the k z = 
plane for a four-band model with p = O.lnm" 3 and 
h = 0.0275Ry. This value of h solves the mean-field equa- 
tion at T = 0, Nmu — lnm -3 , andpA^ = 0.15eVnm 3 , which 
is the upper experimental estimateliS for the exchange cou- 
pling constant. The magnetization orientation is along the 
(100) direction. The minority heavy-hole band is empty at 
this value of h. 

Finally we turn to a series of illustrative calcula- 
tions intended to closely model the ground state of 
Ga.95Mn.05As. For this Mn density and the smaller val- 
ues of Jpd favored by recent estimates, h(T — 0) = 
JpdNMnJ ~ O.OlRy. This value of h is not so much 
smaller than the spin-orbit splitting parameter in GaAs 
(A so — 0.025 Ry), so that accurate calculations require 
a six band model. Even with x fixed, our calculations 
show that the magnetic anisotropy of Ga.95Mn.05As fer- 
romagnets is strongly dependent on both hole density 
and strain. The hole density can be varied by changing 
growth conditions or by adding other dopants to the ma- 
terial, and strain in a Ga.95Mn.05As film can be altered 
by changing substrates as discussed previously. The cu- 
bic anisotropy coefficients (in units of energy per vol- 
ume) for strain-free material are plotted as a function of 
hole density in the inset of Fig. [| the main plot shows 
the coefficients in units of energy per particle. Over the 
density range p < 0.05nm~ 3 , four and six band mod- 
els are in good agreement. We see from this result that 
the asymptotic low density region where the anisotropy 
energy varies as p 5 / 3 holds only for p < 0.005nm~ 3 at 
this value of h. The easy axis is nearly always deter- 
mined by the leading cubic anisotropy coefficient Kf a , 
except near values of p where this coefficient vanishes. 
As a consequence the easy-axis in strain free samples is 
almost always either along one of the cube edge direc- 
tions (Kf a > 0), or along one of the cube diagonal di- 
rections (Kf a < 0). Transitions in which the easy axis 
moves between these two directions occur twice over the 







range of hole densities studied. (Similar transitions occur 
as a function of h, and therefore temperature, for fixed 
hole density.) Near the hole density 0.01 nm~ 3 , both 
anisotropy coefficients vanish and a fine-tuned isotropy is 
achieved. The slopes of the anisotropy coefficient curves 
vary as the number of occupied bands increases from 
1 to 4 with increasing hole density. This behavior is 
clearly seen from the correlation between oscillations of 
the anisotropy coefficients and onsets of higher band oc- 
cupations, plotted in Fig. |Io|. 

Six-band model Fermi surfaces are illustrated in 
Figs. fu||l3| by plotting their intersections with the k z = 
plane at p = 0.1nm~ 3 for the cases of (100), (110), and 
(111) ordered moment orientations. Comparing Fig. Q 
and Fig. which differ only in the band model em- 
ployed, we see that there is a marked difference between 
the majority-spin heavy hole bands in four and six-band 
cases. For the six band model, quasiparticle dispersion is 
particularly slow, leading to large Fermi radii along the 
(110) directions. The large and more anisotropic mass 
is a consequence of mixing with the split-off hole bands. 
As we see from Figs. |l2] and |l3|, this effect occurs for 
all ordered moment orientations, although the details of 
the small minority band Fermi surface projections change 
markedly. The dependence of quasiparticle band struc- 
ture on ordered moment orientation, apparent in compar- 
ing these figures, should lead to large anisotropic mag- 
netoresistance effects in IIIi-^Mn^V ferromagnets. We 
also note that in the case of cube edge orientations, the 
Fermi surfaces of different bands intersect. This prop- 
erty could have important implications for the decay of 
long- wavelength collective modesEf . 
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FIG. 10. Six-band model band densities as a function of 
the total hole density; h — O.OlRy. 



0.10 



0.05 



0.00 



-0.05 




-0.10 



-0.10 



-0.05 0.00 



0.05 0.10 



FIG. 11. Six-band model Fermi surface intersections with 
the k z — plane for p = O.lnm -3 and h = O.OlRy. This figure 
is for magnetization orientation is along the (100) direction. 
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FIG. 9. Six-band model, h — O.OlRy, results for the cu- 
bic magnetic anisotropy coefficients in units of Rydberg per 
particle (main plot) and in kj per m 3 (inset). 
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FIG. 12. Six-band model Fermi surface intersections with 
the k z — plane for the parameters of figure [ll] and magne- 
tization orientation along the (110) direction. 
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FIG. 13. Fermi surface intersections with the k z = plane 
for the parameters of figure [n] and magnetization orientation 
along the (111) direction. 

In Fig. |IJ we present mean-field theory predictions 
for the strain-dependence of the anisotropy energy at 
h = O.OlRy and hole density p = 0.35nm -3 . Accord- 
ing to our calculations, the easy axes in the absence of 
strain are along the cube edges in this case. This cal- 
culation is thus for a hole density approximately three 
times smaller than the Mn density, as indicated by re- 
cent experiments. The relevant value of eo depends on 
the substrate on which the epitaxial Ga.95Mn.05As film 
is grown, as discussed in Section IV. The most important 



conclusion from Fig. [14] is that strains as small as 1% are 
sufficient to completely alter the magnetic anisotropy en- 
ergy landscape. For example for (Ga,Mn)As on GaAs, 
eo = —.0028 at x = 0.05, the anisotropy has a rel- 
atively strong uniaxial contribution even for this rela- 
tively modest compressive strain, which favors in-plane 
moment orientations, in agreement with experiment. A 
relatively small (^1 kJ m -3 ) residual plane- anisotropy 
remains which favors (110) over (100). For x = 0.05 
(Ga,Mn)As on a x — 0.15 (In,Ga)As buffer the strain is 
tensile, eo = 0.0077, and we predict a substantial uni- 
axial contribution to the anisotropy energy which favors 
growth direction orientations, again in agreement with 
experiment. For the tensile case, the anisotropy energy 
changes more dramatically than for compressive strains 
due to the depopulation of higher subbands, as shown in 
Fig. |lj. At large tensile strains, the sign of the anisotropy 
changes emphasizing the subtlety of these effects and the 
latitude which exists for strain-engineering of magnetic 
properties. 
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FIG. 14. Energy differences among (001), (100), (110), and 
(111) magnetization orientations vs. in-plane strain eo at 
h — 0.01 Ry and p = 3.5 nm -3 . For compressive strains 
(eo < 0) the systems has an easy magnetic plane perpendic- 
ular to the growth direction. For tensile strains (eo > 0) the 
anisotropy is easy-axis with the preferred magnetization ori- 
entation along the growth direction. The anisotropy changes 
sign at large tensile strain. 
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FIG. 15. 3rd and 4th band densities for (100), and (111), 
and (001) magnetization orientations vs. in-plane strain eo at 
h = 0.01 Ry and p = 3.5 nm~ 3 . Curves for the (110) magne- 
tization orientation (not shown here) are similar to those for 
the (100) orientation. 



VI. DISCUSSION 

We first comment on the implications of the consider- 
ations described in this paper for the interpretation of 
current experiments. The hysteretic effects which re- 
flect magnetic anisotropy have been studied most ex- 
tensively for the highest T c samples currently avail- 
able. These mean-field-theory predictions depend on 
three phenomenological parameters, Nmu which is sam- 
ple dependent but accurately known, J p d which should 
be nearly universal for a given III-V host compound but 
is less accurately known, and the hole density p which 
is sample dependent and not accurately known. Values 
of Jpd and p must be inferred from experiment, some- 
times by comparison with theoretical pictures which are 
not yet fully developed. The reliability of J p d and p esti- 
mates is improving and presumably will continue to im- 
prove. It now seems clear that the values of J p d and 
p in current high T c samples are such that several hole 
bands are partially occupied in the ferromagnetic ground 



state. In this case, we see from Fig. || that our calcu- 
lations predict cube edge easy axes which include the 
growth direction. The anisotropy energies are typically 
~ 10 -6 Rynm~ 3 ~ lkJm~ 3 . Similar anisotropy ener- 
gies are produced by strains as small as eo ~ 0.001 and 
more typical strains produce larger anisotropy energies. 
An important conclusion from this work is that strain 
contributions to the anisotropy will not normally be neg- 
ligible. 

We believe that the in-plane easy-axis observed in 
Gai^Mn^As films grown on GaAs is a consequence of 
compressive strain in the magnetic film which dominates 
the cubic band anisotropy energy. When Gai-^Mn^As 
is grown on (In,Ga)As, the lattice-matching strain is ten- 
sile, reinforcing the growth direction easy-axis anisotropy 
of strain free samples. (As discussed earlier, the signs of 
both strain and cubic band contributions to strain change 
when the light hole bands are depopulated.) We note 
that in our calculations, the cubic band anisotropy is 
almost always dominated by the fourth cubic harmonic 
coefficient. Given this, it follows that only the cube edge 
easy axes, which includes the growth direction axis, and 
cube diagonal axes, which are not in the film plane, are 
possible without strain. 

Since the local moments are fully polarized in the fer- 
romagnetic ground state, it is easy to estimate the sat- 
uration moment M s ~ Nmti9lHbJ , leading to the rel- 
atively small numerical value [j,qM s ~ 0.05T. It follows 
that the growth direction orientation magnetostatic en- 
ergy ~ ^oMg ~ O.lkJm -3 , considerably smaller than 
the magnetocrystalline anisotropy coefficient. Even when 
several hole bands are occupied with competing spin- 
orbit interactions, these materials have relatively large 
magnetic hardness parameters. Unlike the case of metal- 
lic thin film ferromagnets which have much larger satura- 
tion moments, the magnetostatic shape anisotropy plays 
a minor role in the dependence of total energy on moment 
orientation. The small saturation moment, will also tend 
to lead to large domain sizes and square easy-axis direc- 
tion hysterisis loops, as seen in experiment. 

Coercivities can be estimated^ from the anisotropy 
fields defined by 

K 

fi H a ~ fJ-oM s — j . (37) 

For hard magnetic materials, anisotropy fields are much 
larger than saturation magnetizations. The itinerant 
field places an upper bound on and is expected to scale 
with the coercivity. Our calculations suggest that the 
coercivity in ferromagnetic samples with a single par- 
tially occupied hole band will be immensely larger than 
the coercivity of current samples. Such samples could 
be fabricated, for example, by adding donors such as Si 
to current samples, further compensating the Mn accep- 
tors. According to mean-field-theory, this modification 
in the sample preparation procedure will lower the fer- 
romagnetic critical temperature, and at the same time, 
increase the anisotropy energy. 
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Finally we conclude with a few words of caution. This 
theory of magnetic anisotropy has three principle limita- 
tions: i) it is based on a mean-field theory description of 
the exchange interaction between localized spins and va- 
lence band holes, ii) it neglects hole-hole interactions, and 
iii) it doesn't account for disorder scattering of the itin- 
erant holes. Confronting these weaknesses would in each 
case considerably complicate the theory and we feel it 
is appropriate to seek progress by comparing the present 
relatively simple theory with experiment. Nevertheless it 
is worthwhile to speculate on where and how the theory 
may be expected to fail. 

Mean-field theory should be reliable when the range of 
the hole mediated interaction between localized spins is 
long and the spin-stiffness parameter that characterizes 
the energy of long-wavelength spin-fluctuations is_suffi- 
ciently large. Considerations of this type suggested that 
mean-field theory will fail at high temperatures unless 
the ratio of the hole density to the localized spin density 
is small and the field h is not too large compared to the 
Fermi energy. Since p is typically smaller than Nmu be- 
cause of anti-site defects in low-temperature MBE growth 
samples, mean-field theory is likely to be reasonable at 
least at T — in many (III,Mn)V ferromagnets. 

Hole-hole interactions will clearly tend to favor the fer- 
romagnetic state by countering the band-energy cost of 
the spin polarization. Because of strong spin-orbit cou- 
pling in the valence band, estimates based on many-body 
calculations for electron gas systems may be of little use 
in estimating the importance of this effect more quan- 
titatively. Work is currently in. progress which should 
shed more light on this mattertZl. Nevertheless, it seems 
likely that these interactions will not have an overriding 
importance, at least in large hole density samples. 

Finally we come to disorder. It is clear from experi- 
ment, that disorder does not have a qualitative impact 
on free-carrier mediated ferromagnetism even when those 
free carriers have been localized by a random disorder 
potential. It seems likely that disorder will destroy the 
ferromagnetic state, only when the localization length 
becomes comparable to the distance between localized 
spins. On the other hand, since elastic disorder scatter- 
ing will mix band states with different orientations on 
the Fermi surface, it also seems clear that a reduction in 
magnetic anisotropy energy must result. Indeed the co- 
ercivities that follow from our anisotropy energy results 
appear to be larger than what is observed. As far as we 
are aware, no theory of this effect exists at present. 
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APPENDIX 

In the literature, different representations are used for 
the four band and six band model Kohn-Luttinger Hamil- 
tonians. In the interest of completeness and clarity, this 
appendix specifies the expressions on which our detailed 
calculations are based. Detailed derivations of k • p per- 
turbation .theory for cubic semiconductors can be found 
elsewhere. E3 

The k = states at the top valence band have p-like 
character and can be represented by the I = 1 orbital 
angular momentum eigenstates \mi). In the coordinate 
representation we can write 

(r|mj = 1) = --j=f{r){x + iy) 



(r\mi = 0) = f(r)z 
(r\mi = -1) = —f{ r ){x - iy). 



(38) 



The Kohn-Luttinger Hamiltonian for systems with no 
spin-orbit coupling, Hz,, is written in the representation 
of the following combinations of |mj) 



1 



\X) = -j={\m l = 

\y) = ^{\m l = 

\Z) = \mi = 0) . 



-i) - \m = i» 
-1) + \mi = 1» 



It reads 



/ Akl + B{kl + kl) 



Ck x ky 



Cky k x 
Ck z k x 



Akl + B(kl + kl) 



Ck z ky 



Akl 



(39) 



Ck x k z 
C^jky k z 
f B(kl - 

(40) 



where 
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B 



C 



-^(71 -272), 

73, 



(41) 



m is the bare electron mass, and 71, 72, and 73 are 
the phenomenological Luttinger parameters. To include 
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spin-orbit coupling we use the basis formed by total an- 
gular momentum eigenstates \j,rrij): 

\l) = \j = 3/2,m s =3/2) 

|2) = 1^=3/2,^- = -1/2) 

|3) = 1^ = 3/2^ = 1/2) 

|4) = |i = 3/2,m, = -3/2) 

|5> = |i = l/2,m, = l/2) 

\6) = \j = l/2,m s = -l/2) (42) 

The basis (^2|) is related to the orbital angular momen- 
tum (mi = 1, 0, —1) and spin (a =|, {) eigenstates by 



|l) = |m, = l,T) 

12 



i|m, = -l,t) + ^||m I = 0,i> 

|3) = -^K = i,|) + yf K = o,T) 
|4) = K = -ia) 

|5> = -^|m/ = 0,T> + Y / fl^ = 14> 

|6) = iK = o,i)-yf N = -i,T) 



(43) 



.-^(|x,t) + i|r,t)) 



|2> = 



7! (|X ' T) 



^4) 



|3) = --^(|x,t) + i|r,t)) + y||^T) 

| 4 ) = i=(|x ) 4)-i|F,|)) 

|5) = ~(|x,i) + »|y > i))--^|z > T) 

\6)=-^(\X,V-i\Y,V) + -^\Z,l) 



(44) 



The six band model Kohn-Luttinger Hamiltonian, Hi,, 
in the representation of vectors ([42]) is 



H L = 



( Hhh -c 

-c* Hi h 

-b* 

b* 



b* 
V2 



\c*V2 



b\/3 
' V2 

-d* 



-b 








6 


7~Lih 


— c 


— c* 


7~lhh 


d* 


-cV2 


6*\/3 


6 


%/2 





%/2 

fc*V3 



cV2 \ 
-d 

%/2 





(45) 



In the matrix ( f45| ) we highlighted the four band model 
Hamiltonian. The Kohn-Luttinger eigenenergies are 
measured down from the top of the valence band, i.e., 
they are hole energies and we use the following notation: 



n hh = — [(71 + i2){k 2 x + k 2 y ) + ( 7l - 2 l2 )k 2 z 

Hih = [(71 - 70(*£ + k 2 y ) + (71 + 2 72 )fc z 2 
Im u 

h 2 

Hso = 7^1l{ k l +K + k z) + A *° 
2m y 

b = ^^—J3k z (k x - iky) 
m 

c = [l2(kl - k 2 y ) - 2ij 3 k x k y ] 



2m 
V2h 2 
2m 



l2 [2kl-{kl + kl)] 



(46) 



The four band Kohn-Luttinger Hamiltonian can be di- 
agonalized analytically and yields a pair of Kramers dou- 
blets with eigenenergies 



In the spherical approximationta (72,73 — > 7 = O.672 + 
O.473), the top of the valence band consists of two doubly 
degenerate parabolic bands with effective masses m^ — 
771/(71 - 27) ~ 0.498m and m, = 777/(71 + 27) - 0.086m. 
(An approximation in which the light hole bands, which 
have a much smaller density of states, are ignored is ad- 
equate for some purposes.) 

The four band and six band representations for the 
hole spin-operator components read 









1 

2V3 





1 

%/6 




\ 








1 

3 


1 

2V3 




1 

3V2 





1 

2^3 


1 
3 













1 

3V2 





1 













1 

V6 


1 


1 

3%/2 













1 

6 








1 

3V2 


1 

Ve 


1 
() 




/ 


( 





1_ 

2v^3 







1 

Ve 











1 

3 


1 

2\/3 




1 

3V2 





1 

2V3 


1 

3 













1_ 

3V2 





1 

2V3 













1 

71 


1 

75 


1 

3V2 













1 

6 


V 





1 

3V2 


1 




1 

6 
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/ 1 

2 











__L 

2 





u 3 





u 3 






V2 
3 



v2 

3 





_L 

6 



(47) 
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Host 


7 < 100 > 


7 < 110 > 


7 < 111 > 




1? 


GaAs 


5.965 


5.088 


4.639 


-3.509 


-4.24 


InAs 


13.207 


10.854 


9.705 


-9.412 


-9.84 



TABLE I. High symmetry direction moment-orientation 
dependent average Luttinger parameters and their cubic har- 
monic expansions for GaAs and InAs based IIIi-^Mn^V fer- 
romagnetic semiconductors. These parameters specify the 
T — magnetic anisotropy energy in the limit of large 
spin-orbit splitting and large exchange coupling parameter 
( J p d) or small hole density p. 
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